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ABSTRACT 

SCUBA-2 is the next-generation replacement for SCUBA (Sub-millimetre Common User Bolometer Array) on 
the James Clerk Maxwell Telescope. Operating at 450 and 850 /im, SCUBA-2 fills the focal plane of the telescope 
with fully-sampled, monolithic bolometer arrays. Each SCUBA-2 pixel uses a quarter-wave slab of silicon with 
an implanted resistive layer and backshort as an absorber and a superconducting transition edge sensor as 
a thermometer. In order to verify and optimize the pixel design, we have investigated the electromagnetic 
behaviour of the detectors, using both a simple transmission-line model and Ansoft HFSS™, a finite-element 
electromagnetic simulator. We used the transmission line model to fit transmission measurements of doped 
wafers and determined the correct implant dose for the absorbing layer. The more detailed HFSS modelling 
yielded some unexpected results which led us to modify the pixel design. We also verified that the detectors 
suffered little loss of sensitivity for off-axis angles up to about 30°. 
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1. INTRODUCTION 

In SCUBA-2 1 two detector arrays will observe simultaneously at wavelengths (A) of 850 and 450 /im. Each array 
consists of four 40 x 32 sub-arrays of square pixels. Each pixel comprises a silicon slab which has a metal layer, 
or backshort, behind it. For the interface between free space and the front of the slab, the reflection coefficient is 
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where n is the refractive index of the slab. We assume n = 3.388 for cold silicon, giving R — 30%. However, we 
can reduce the reflectivity (and thus maximize the detector efficiency) if the thickness of the slab is chosen so 
that the radiation reflected by the backshort at the back of the slab interferes destructively with the radiation 
reflected by the top of the slab. In both cases, the radiation undergoes a 180° phase change on reflection, so 
the slab thickness is chosen to be an odd multiple of A/4 in the slab material. This arrangement is known as a 
Salisbury screen and is used in aviation for reducing radar reflections. The operation of the detector is shown 
in Figure 1. At sub- millimetre wavelengths and low temperatures, silicon has negligible dissipation. Thus, the 
slab has an implanted resistive layer at the top to dissipate the energy of the radiation, heating up the pixel. 
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To maximize the detector efficiency, we need to choose the implantation dose so that the surface resistance is 
R s ~ 377 fi/b, matching the wave impedance of free space (Zq = 120 tt O). The temperature rise due to the 
absorption of radiation is detected by a transition-edge sensor (TES) on the back side of the quarter-wave slab. 
It is this TES that acts as the reflective backshort. 

The pixel geometry is repeated across the detector wafer to form one of SCUBA-2's 40 x 32 sub-arrays. Each 
sub-array is bump-bonded to a wafer that carries the SQUID multiplexing (MUX) circuitry. The indium bump 
bonds provide electrical and thermal connections, as well as mechanical strength. 

The size of an individual quarter-wave slab (1055 fira, allowing for the wall thickness and the gap between the 
walls and slab) is on the order of the wavelength of the incident radiation (850 and 450 /iin). This means that to 
design the pixels for maximum efficiency we need a model of the electromagnetic behavior of the detector. Wc 
need answers to the following questions. How should we design the pixels to maximize the detector efficiency, and 
how much variation in the detector parameters can we tolerate? How does the efficiency of the detector depend 
on the wavelength, polarization, and incident angle? To answer these questions, we carried out electromagnetic 
modeling. We approached the problem in two ways. First, we modeled the detectors as a simple transmission 
line. This model is equivalent to an infinite dielectric slab absorber and does not take into account any effects 
due to the finite size of the pixels, the supporting grid, or the periodicity of the array. Then, we used a com- 
mercial software package, Ansoft HFSS™ (High-Frequcncy Structure Simulator; see http://www.ansoft.com), 
to model the behaviour of the detectors more realistically. HFSS is a three-dimensional finite-clement method 
electromagnetic simulator. 
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Figure 1. Operation of the SCUBA-2 pixels. The detectors consist of a metal reflector and a resistive layer spaced a 
quarter wavelength apart by a slab of silicon. This arrangement is known as a Salisbury screen. The emergent waves (ei) 
interfere destructively with the reflected wave (R) and constructively with the incident wave. This maximizes the energy 
dissipated in the resistive layer at the top of the quarter-wave slab and minimizes the amount reflected. 



2. TRANSMISSION LINE MODEL 

For a simple model we can represent the vacuum, silicon, metal sandwich as a transmission line with corresponding 
impedances (see Figure 2). Because we are assuming that the pixels extend laterally to infinity, we neglect any 
effects due to the finite size of the pixels or the periodicity of the array. These effects will be addressed in 
Section 3. Using this simple transmission line model, the voltage transmission coefficient is 
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Figure 2. SCUBA-2 quarter-wave slab represented by a simple transmission line model. The slab has a thickness d. The 
surface impedances of its front and rear faces are Z s and Zt- The dielectric of refractive index n is represented by an 
impedance Zsi = Zo/n, where Zo is the impedance of free space. For the SCUBA-2 pixels we have R s = 377 Q/o and 
Z t -> 0. 
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Figure 3. Dependence of quarter-wave slab absorption coefficient on slab thickness predicted from the simple transmission 
line model for wavelengths of 850 (left) and 450 /J,m (right). The quarter- wave slab has a resistive layer on the front with 
R s = Z . 



and the voltage reflection coefficient is 
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where z s = Z s j 'Zq and z t = Z t / Z$. The corresponding power transmission, reflection, and absorption coefficients 
are T = |t| 2 , R = |r| 2 , and A = 1 — T — R. respectively. In the general case, these expressions are complicated. 

To model the behaviour of the SCUBA-2 detectors, we let z t — * in Equation 2. This means that we are 
assuming that the back of the slab is covered by a perfect conductor (the backshort). Figures 3 and 4 show how 
the predicted absorption coefficients depend on the slab thickness and the wavelength of the incident radiation. 
We are assuming that the implanted absorbing layer is purely resistive. Note that the |A slab of the 450 /zm 
detectors has greater wavelength selectivity than the |A slab of the 850 /jm detectors. 

2.1. Angular Dependence 

In the current design, the SCUBA-2 optics are not telecentric. Pixels at the edges of the arrays will have rays 
incident on them at angles up to about 18°. Thus, we must ensure that the efficiency of the detectors does not fall 
off too rapidly with incident angle. Otherwise, pixels near the edges of the arrays will be insensitive. We used the 
simple transmission line model to investigate the angular dependence of the detector efficiency. We generalized 
the model to oblique incidence by making the substitutions shown in Table l. 2 ' 3 For normal incidence, the 
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Figure 4. Wavelength dependence of absorption coefficient of quarter-wave slab predicted from the simple transmission 
line model. Left: Slab thickness d = 62.7 /im (850 detector). Right: Slab thickness d = 99.6 fim (450 /im detector). 
The quarter- wave slab has a resistive layer on the front with R 3 = Zq. 



Table 1. Substitutions to generalize transmission line model to case of oblique incidence. ^t(Zsi) 
is the real part of Zsi- 
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two linear polarizations are equivalent. However, for oblique incidence, we must treat the s-polarization (E-field 
parallel to the plane of incidence) and the p-polarization (E-field normal to the plane of incidence) separately. 

Using the substitutions in Table 1 with the transmission line model we calculated the angular dependence of 
the efficiencies of the SCUBA-2 detectors for s- and p-polarizcd radiation. The results are shown in Figure 5. 
We found that the |A slab of the 450 fim detectors has greater polarization-dependence than the jX slab of the 
850 fim detectors. 
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Figure 5. Angular dependence of absorption efficiency of SCUBA-2 850 fim (left) and 450 fim (right) detectors for s- 
and p-polarized radiation, predicted from the simple transmission line model. 
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Figure 6. Angular dependence of absorption efficiency of SCUBA-2 850 /J,m (left) and 450 fim (right) detectors for s- and 
p-polarized radiation, predicted from the simple transmission line model assuming a 10% frequency bandpass (modeled 
as a double Gaussian). The on-axis efficiency is degraded to 99.14% of its monochromatic value for the 850 /im detectors 
and 93.35% for the 450 /im detectors. 



2.2. Effect of Filter Bandpass 

The radiation incident on the SCUBA-2 detectors will not be monochromatic. To investigate the effect of the 
finite bandpass of the SCUBA-2 filters, we smeared the results of the transmission line model over a double- 
Gaussian filter function with a frequency bandpass of 10%. The effect of smearing over the filter bandpass is 
slight in the case of the 850 /im detectors. However, for the 450 /xm detectors the on-axis efficiency is degraded 
to 93.35% of its monochromatic value. The results are shown in Figure 6. The difference between the 850 and 
450 /im detectors is not unexpected. As was shown in Figure 4, the 450 /*m detectors have greater wavelength 
selectivity. 

2.3. Implant Dose 

Our first application of the transmission line model to real data was to determine the correct implant dose for 
the implanted layer that absorbs the incident radiation. Wafers were implanted on the front side with different 
doses. The implantation process is described elsewhere in these proceedings. 4 The submillimetre transmission 
of the wafers was measured at 1.5 K as a function of wavenumber (a) using a Fourier transform spectrometer 
at Cardiff University. To model the transmission of these implanted wafers, we let Zt — ► oo in Equation 2. We 
found that the model gave a satisfactory fit only if we added another component to the model — an overall 
normalization factor — and allowed this to vary freely. The fits are consistent with the implanted layer being 
purely resistive rather than metallic. Assuming we have a purely resistive implanted layer, as long as the correct 
thickness is chosen for the quarter-wave spacer the absorption efficiency does not depend strongly on the surface 
resistance R s (see Figure 7). 

3. FINITE-ELEMENT MODELING WITH ANSOFT HFSS 
3.1. Description of the Software 

The simulations described here were carried out using HFSS version 8.5. HFSS includes a three-dimensional 
drawing program that can be used to construct a geometrical model. In other parts of the software, the user 
specifies the material properties, boundary conditions, and electromagnetic excitation of the model. HFSS 
divides the volume of the model into a mesh of tetrahedra and solves for the electric field on the mesh in the 
frequency domain. The problem is converted into a set of linear equations, represented by a matrix. A field 
solution is obtained by inverting the matrix for the initial mesh. The mesh is then refined, based on the local 
field strength. The field solution from the refined mesh is compared with the previous solution. This process of 
refinement and solving continues until the user-specified convergence criteria are met. The convergence criteria 
depend on the type of excitation. For the models described here the model is excited with an incident plancwave 
and convergence is evaluated from changes in the total scattered energy, calculated from the E-field. 
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Figure 7. Dependence of absorption efficiency of SCUBA-2 850 (left) and 450 /im (right) detectors on surface resistance 
R s predicted from the simple transmission line model. 

HFSS includes tools for visualizing the field solutions and calculating physical quantities from them. These 
may be run from scripts to calculate quantities such as the reflection and absorption coefficients automatically. 
We used the Ansoft Optimetrics™ package to run HFSS repeatedly on models while varying geometric or material 
parameters. 

3.2. Description of the Models 

3.2.1. Method 

We simulated an infinite array of pixels by modeling a single square pixel and applying periodic boundary 
conditions to opposite walls. The models are excited by an incident plane wave. HFSS solves for the E-ficld 
in the model. From the E-field we can calculate the physical quantities that interest us. These are the power 
absorbed in the resistive layer, the total incident and reflected power, and the power flowing in or out of the unit 
cell through the side walls. The power absorbed in the implanted resistive layer is calculated from the field on 
the surface: 

P abs = ^-J s (Exh) 2 dS (4) 

where R s is the surface resistance of the implanted layer and n is the unit normal to the surface. We have 
assumed R s = 377 fi/n for all of the models. The other powers are found by integrating the Poynting vector 
over the appropriate surfaces. 

Our basic pixel model consists of an approximation to a unit cell of the array (see Figure 8). Vertical 
boundaries run through the centres of the silicon walls. The walls are 50 /im thick, so that the unit cell contains 
the half-width of 25 /im. Periodic boundary conditions are applied to opposite side walls. In the original pixel 
design there was a supporting grid 381 /im high on top of the walls. This was later removed from the design (see 
Section 4.3). The results presented here are for models without this grid, unless otherwise stated. The volume 
of the model extends lj wavelengths above the top of the walls (or grid, if present). Inside the walls is the 
quarter-wave silicon slab with a resistive layer on its top surface. Under the slab is a perfectly conducting layer, 
to simulate the superconducting TES. 

3.2.2. Selection of boundary conditions 

For the top of the model where the waves come in, we have a choice of radiation boundary or perfectly matched 
layer (PML). Both of these are designed to emulate free space by absorbing incident radiation. The radiation 
boundary will absorb normally-incident radiation but it has a finite reflectivity for large angles of incidence. The 
PML is a lossy, anisotropic dielectric that will also absorb radiation incident at large angles. Because we wanted 
to illuminate the pixels with radiation incident at oblique angles, we chose a PML for the top of the model. 

On the side walls, we apply periodic boundary conditions for consistency with the plane wave excitation. The 
periodic boundary conditions on the pixel walls result in an asymmetric matrix, which increases the problem 
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Figure 8. Pixel model used in HFSS simulations. 



size. However, neither radiation boundaries or PMLs arc suitable for the side walls and yield results inconsistent 
with energy conservation. 

3.2.3. Error Analysis 

Because HFSS provides no direct measure of the errors in the results, we needed to find some way to check the 
quality of solutions. The results can be checked for consistency by accounting for all of the incident power flowing 
into the unit cell. The integrated flux flowing through all the walls of the model plus the power absorbed in the 
resistive layer of the quarter-wave slab should add up to the incident power. Any difference should be small for a 
high-quality solution. In the plots below, the error bars on the absorption and reflection coefficients are derived 
from this discrepancy in the total power. While they are not useful for deriving statistical uncertainties, they 
are shown so that the relative accuracy of the solutions can be compared. Another way to evaluate the solution 
is to calculate the incident power flowing through the upper boundary and compare it with the Poynting flux 
integrated over the surface. The results should be slightly different due to the finite size of the mesh elements. 
In HFSS, for plane wave excitation, the default setting is to apply a free-space plane wave with a peak electric 
field of 1 V m _1 . From this we can calculate the incident flux as 

P mc = [ S.dA = -ExHi cos # = - < f^-E 2 A cos 8 (5) 

where A is the area of the top surface of the model, and 9 is the angle of incidence. The factor of 1/2 arises from 
the fact that the HFSS field calculator uses peak phasors to represent the fields. We found that the incident flux 
calculated from Equation 5 sometimes differed from the integrated Poynting flux by up to 10%. Thus, in the 
results given below, the incident flux is derived from Equation 5. 
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Dependence of Absorption Coefficient on Incident Angle 
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Figure 9. Comparison of results from HFSS (points) and the simple transmission line model (lines) for a Salisbury screen. 



In order to check the accuracy of the HFSS solutions we simulated the 850 and 450 /im detectors as infinite 
slabs and compared the results with those from the transmission line model. Figure 9 shows good agreement for 
both, although the 450 (im efficiency is underestimated and has greater uncertainty. This is because the shorter 
wavelength requires a finer mesh for the same precision. The ultimate limit on the precision of the field solution 
is memory. This limits the total number of tetrahedra in the mesh. The version of HFSS we used is a 32-bit 
program. This means that it will fail if the mesh grows to more than 4 GB. In order to obtain the best precision 
possible, we pushed the mesh refinement until this limit was reached when solving our models. 

4. RESULTS OF THE HFSS MODELING 

4.1. Detection Efficiency 

We used the Ansoft Optimetrics package to investigate how the absorption coefficient varied with the thickness 
of the quarter-wave slab. The results are shown in Figure 10. For the 850 /im detectors the absorption coefficient 
is greatest with a slab thickness of ~ 65 /im. From the assumed value of the refractive index (n = 3.388) we 
would expect the maximum to occur at a slab thickness of 62.7 /*m which corresponds to A/4 in silicon. For 
the 450 /im detectors the maximum occurs close to the expected value of 100 /jm (3A/4). In each case the peak 
monochromatic detection efficiency was 94%. 

4.2. Angular dependence 

We investigated the off-axis response of both the 850 and 450 /im detectors with the Ansoft Optimetrics package. 
Neither detector suffers much loss in sensitivity for incident angles up to 30° (see Figure 11). However, this 
sensitivity to oblique radiation does make the detectors susceptible to stray light. 
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Figure 10. Predicted variation of detector efficiency with quarter-wave slab thickness for 850 /iiri (left) and 450 /im 
(right) detectors from HFSS simulations. The vertical dashed lines correspond to the expected thicknesses for maximum 
efficiency. 
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Dependence of Detector Efficiency on Incident Angle 
for 450 /mi Detectors with s-Polarized Incident Wave 
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Figure 11. Predicted variation of 850 /im (top) and 450 /mi (bottom) detector efficiency with angle of incidence from 
HFSS simulations. Left: s-polarization. Right: p-polarization. 



4.3. Effect of Supporting Grid 

The original design for the SCUBA-2 pixels had a silicon grid supporting the detector wafers. A model that 
includes this grid is shown in Figure 12. This grid was essentially an extension of the side walls by 381 /im, a 
standard wafer thickness. The purpose of the grid was to provide mechanical rigidity to the detector array. Since 
silicon is essentially transparent at these wavelengths, its effect on the detector performance was thought to be 
insignificant. However, the HFSS modeling revealed that the grid acted as a dielectric waveguide, both reducing 
the detector efficiency and changing the optimum slab thickness (see Figure 13). 
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Figure 12. Pixel model corresponding to original pixel design. Note that the silicon support grid has been removed from 
the design and is omitted from simulations described here, unless otherwise stated. 

The mechanical support grid also affects the angular dependence of the detector efficiency (see Figure 14). 
The grid structure behaves like a two-dimensional diffraction grating when the wavelength is smaller than the 
grid spacing. At both wavelengths (450 and 850 fxva) the array allows diffraction orders to propagate. This 
means that the non-absorbed reflected field is a superposition of the zeroth order (specular) component plus all 
the diffraction orders that can be generated. These depend on the frequency, angle of incidence and the grid 
spacing (i.e. the array pitch). HFSS models this behaviour in a region that is close to the array, i.e. in its 
near field. This means that we must be careful when comparing the simulations with experimental data. When 
we measure the reflected flux wc can be in a far-field configuration where only the zeroth order is detected. In 
this case the HFSS results have to be processed to calculate the far-field in the periodic boundary case. The 
important point is that the method we have described to calculate the absorption remains correct. 




Figure 13. Predicted variation of detector efficiency with quarter-wave slab thickness for 850 /jm (left) and 450 /an (right) 
detectors from HFSS simulations that include a 381 /im mechanical support grid. The vertical dashed lines correspond 
to the expected thicknesses for maximum efficiency. 
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Figure 14. Predicted variation of 850 /mi (top) and 450 /im detector efficiency with angle of incidence from HFSS 
simulations that include a 381 /im mechanical support grid. The vertical dashed lines represent the diffraction minima 
and maxima for the array pitch (1135 /tm). Left: s-polarization. Right: p-polarization. 



Because of its effect on the detector efficiency and because changes in the processing had made it redundant 
anyway we removed the mechanical support grid from the design. 

4.4. Flatness of quarter-wave slab 

We used the Ansoft Optimetrics package to investigate the effect of the top of the quarter-wave slab not being 
fiat. To the top of the slab we added a wedge, and adjusted the thickness of the slab and that of the wedge, 



keeping the mean thickness of the combined structure constant. This was done so that we could see the effects 
of having the top of the slab not parallel to the backshort, independent of changes in the thickness of the slab. 

The top of the quarter- wave slab can make an angle of up to 1.6° with the backshort before reducing the 
absorption coefficient for the p-polarization by 5%. This corresponds to a deviation in the slab thickness of up 
to ±15 [im across a pixel. 

4.5. Effect of oxide layer 

The fabrication process requires a layer of silicon oxide on the slab and walls as an etch stop. We ran simulations 
with a 1 /jm thick layer of silicon dioxide on top of the walls and the slab. It had no significant effect on our 
results. 

5. CONCLUSIONS 

• The optimum quarter- wave slab thicknesses to attain at least 95% of the maximum efficiency are 66 ±6 ^tm 
and 99 ± 3 /im for the 850 and 450 /j,m detectors, respectively. 

• The optimum surface impedance is 377 fl/n. The absorption is not highly sensitive to the surface 
impedance. In order to attain at least 99% of the optimum efficiency the surface impedance must be 
between 309 and 460 fl/n for the 850 /im detectors and between 308 and 461 il/n for the 450 /xm detec- 
tors. 

• The monochromatic efficiency of the detectors at the optimum slab thickness is 94%. 

• When we take into account the finite bandwidth of the filters, the efficiency of the detectors at the optimum 
slab thickness is 93% at 850 /mi and 88% at 450 /im. 

• The detector efficiency is not very dependent on incident angle up to about 30°. There is no significant 
loss of sensitivity over the range of incident angles in the field of view. 

• The flatness requirement for a decrease in efficiency of less than 5% is that the slab thickness varies across 
a pixel by no more than 15 /im. 

• The presence of a 1 /zm oxide layer on the walls and slab has a negligible effect on efficiency. 

• The removal of the mechanical support grid improves the detector efficiency significantly. 

We conclude from this modeling that the SCUBA-2 detectors will meet the requirement of an absorption 
efficiency of 75% or greater. 
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